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We provide a non-perturbative expression for the hadron production in electron-positron anni- 
hilation at zero temperature in a strongly coupled, large-7V c SU(N C ) field theory with N t <C N c 
quark flavors. The resulting expressions are valid to leading order in the electromagnetic coupling 
constant but non-perturbatively in the SU(N C ) interactions and the mass of the quark. We obtain 
this quantity by computing the imaginary part of the hadronic vacuum polarization function Id, 
using holographic techniques, providing an alternative to the known method that uses the spectrum 
of infinitely stable mesons determined by the normalizable modes of the appropriated fields in the 
bulk. Our result exhibits a structure of poles localized at specific real values of q 2 , which coincide 
with the ones found using the normalizable modes, and extends it offering the unique analytic con- 
tinuation of this distribution to a function defined for values of q 2 over the complex plane. This 
analytic continuation permits to include a finite decay width for the mesons. By comparison with 
experimental data we find qualitatively good agreement on the shape of the first pole, when using 
the p meson parameters and choosing a proper normalization factor. We then estimate the contri- 
bution to the anomalous magnetic moment of the muon finding an agreement within 25%, for this 
choice of parameters. 



I. INTRODUCTION 



Observables like the electromagnetic coupling constant, a(s), and the anomalous magnetic moment of the muon 
have been accurately measured 0, HJ , and have proven to be very sensitive to the effects of Quantum Chromody- 
namics (QCD), which enter through its contribution to the vacuum polarization 0, 0|. A standard quantity used to 
parametrize these QCD effects is the ratio between the following cross sections, 

a(e + e~ — > hadrons) 

R( s ) = , + ~ ; +— s (!) 

cr(e + e — > j 

which is related to the vacuum polarization II(s) by 

R{s) = 127rlmn,(s + ie). (2) 

At large momenta perturbative QCD is justified, and the corrections to R(s) up to second order in the strong coupling 
constant a s have been computed [6]. On the other hand, in the low energy regime, QCD becomes a non-perturbative 
theory and no ab initio QCD calculations for R(s) exist. Therefore, the experimental information on R(s) in this 
regime is plugged to built the theoretical predictions for a(s) and the anomalous magnetic moment of the muon. This 
is crucial in the case of the muon anomalous magnetic moment, where the energy dependence of the electromagnetic 
kernel enhances the low energy QCD contribution. In fact, the theoretical prediction uncertainty is dominated by 
this so-called hadronic contribution and, at the moment, a discrepancy with the experimental measurement persist 
0. Therefore, at low energies, a non-perturbative calculation on a s , accounting for the full dependence of R(s) on 
the quark mass, is mandatory. 

It is in this context that the gauge/gravity correspondence turns out to be a good alternative as a computational 
method. The dual theory of QCD is still unknown, but we can do calculations in SU(N C ) super Yang-Mills (SYM) 
theories coupled to N t flavors of fundamental matter whose bare mass is an adjustable parameter. These theories 
are very different from QCD, but the contribution to R(s) we want to investigate comes from the strongly coupled 
part of QCD. Although the degrees of freedom differ from those of the theory we are working with, this is also a 
strongly coupled theory, and one may hope that certain properties of QCD in this regime could be observed in this 
setting. This supposition has already been used in other contexts to model hadronization [8|, to find the thermal 
hadron spectrum Jj| and to study extensively the strongly coupled quark-gluon plasma (lfj |. 

In this work we compute the hadronic vacuum polarization in this framework at zero temperature, compare its 
behavior with experimental data, and give an estimate of its contribution to the anomalous magnetic moment of the 
muon. Our results are valid to leading order in the electromagnetic coupling constant, but are non-perturbative in 
a s and the mass of the quark. 
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II. R(s) IN THE GAUGE THEORY 

We will obtain R(s) from its expression in terms of the hadronic part of the vacuum polarization function Tlq(q 2 ), 
which at zero temperature is related to its tensorial form H^"(q) by 



Ml 2 ) = 1 , , (3) 
q z d 



where fi and v run from to d, and d + 1 is the number of dimensions of the gauge theory we are working in. 

The matter content of the theories we are working with consists of fundamental and adjoint fields. The former 
include N f flavors of fundamental fermions 1 $ a and scalars $ a , a = 1, . . . , N t , to which we will refer indistinctly as 
'quarks'. To be as close as possible to QCD, we will couple only the fundamental fields to electromagnetism and 
set the electric charge of the adjoint fields to zero, so they do not affect our results for R(s). For this we introduce 
an additional, dynamical, Abelian gauge field (the photon) that couples with strength (electric charge) e to the 
fundamental fields. This amounts to replacing the S'£/(./V c )-covariant derivative D M by T>^ = — ieA^ when acting 
on the fundamental fields, and adding a kinetic term for the photon. In this way we obtain a SU(N C ) x C/(1) EM gauge 
theory whose Lagrangian is 

C = £su(N c ) ~ \l% + eA^J™ , (4) 
with T^v = d^Au — d l jA lll and the electromagnetic current given by 

= * 7ai v& + U>*v^ - % - {v^y $ , (5) 

where a sum over flavor and color indices is implicit. 

In this theory we compute Uf£ v (q) using the following expression in terms of the correlator of two electromagnetic 
currents of quarks, 

U^(q)=i I d d+1 xe~^Q(x a )([J^(x),jr(0)])- (6) 

Thus in order to study R(s) we must in principle calculate the two-point function in (jB]) from the SU(N C ) x [7(1) EM 
theory. However, as first noticed in the terms in the electromagnetic current ([5]) proportional to the photon 
field (implicit in the covariant derivative) lead to higher-order contributions in e for the correlator in ([6]), and can 
thus be ignored to leading order in e. Furthermore we observe that the two-point function of the remaining terms in 
the current can be calculated in the SU (A c ) theory, since again the effects of including the dynamical photon are of 
higher order in e. We therefore conclude that, to leading order in the electromagnetic coupling constant e, R{s) in an 
SU(N C ) x U(1) EM theory is completely determined by the two-point function of the electromagnetic current ([5]) in the 
SU(N C ) theory. This is the key for us to be able to perform this calculation using the gravitational dual description, 
since the dual of the SU(N C ) x U(1) EM is unknown. 



III. HOLOGRAPHIC DESCRIPTION 



We calculate now this correlator in a JV = 2, SU(N C ) SYM theory coupled to N s flavours of fundamental matter. 
To this end, we use the gravitational dual of this theory [l2j|, given by a stack of N f D7 probe branes in the AdS$ x S 5 
background generated by N c D3 branes in the limit N c >> N t . For this system we obtain an analytic result for R(s), 
while for the theory dual to a D4/D6 configuration we could obtain only numerical results, which we do not discuss 
here since they are qualitatively identical to those of the D3/D7 system. 

In the decoupling limit, the D3-brane solution in the string frame takes the form 




where x = (a; 1 , x 2 , x 3 ) and R is the AdS$ radius given by R 4 = 47r(? s jV c ^, with l s and g s the string length and 
coupling constant, respectively. 

According to the gauge/gravity correspondence, string theory on the background above is dual to (3+l)-dimensional, 
maximally supersymmetric Yang-Mills theory at zero temperature. Consider now N t <C N c coincident D7-brane probes 
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that share the x M = (x°, x 3 ) directions with the D3-branes, extend along the radial direction and wrap an S 3 inside 
the S 5 of ([7]). If we place the D7-branes at a fixed distance L from the D3-branes in the plane perpendicular to 
both of them, the configuration will be supersymmetric and this will ensure stability. Under these conditions the 
Ramond-Ramond field sourced by the D3-branes does not couple to the D7-branes. 

In the gauge theory, the D7-branes correspond to introducing N s flavors of fundamental matter [12j. These funda- 
mental fields propagate along a (3 + l)-dimensional defect and have the same mass m q dictated by the distance L 
according to 

The U(Nf) gauge symmetry on the D7-branes is a global flavor symmetry of the gauge theory. If this theory 
is coupled to electromagnetism as explained above, then the electromagnetic current J™ is the conserved current 
associated to the overall U(l) C C/(iV f ). At strong 't Hooft coupling and large N c , each conserved current of the gauge 
theory is dual to a gauge field on the gravity side [Hjj . To describe the gauge field dual to J™, let A m , m — 0, . . . , 7, 
be the gauge field associated to the overall U(l) gauge symmetry on the D7-branes. Upon dimensional reduction on 
the 3-sphere wrapped by the D7-branes, A m gives rise to a massless gauge field {A^, A r }, three massless scalars, and 
a tower of massive Kaluza-Klein (KK) modes. All these fields propagate on the five non-compact dimensions of the 
D7-branes. Following [l4||, we will work in the gauge A r = 0, and furthermore we will set to zero the scalars and the 
higher KK modes, since they are not of interest here, and unlike in the non-zero temperature case (Isj . this can be 
done consistently. The gauge field A^ is the desired dual to the conserved electromagnetic current J™ of the gauge 
theory. 

According to the gauge/gravity correspondence, correlation functions of J™ can be calculated by varying the string 
partition function with respect to the value of A^ at the boundary r — > oo of the spacetime |13j . Under the present 
circumstances, the string partition function reduces to e tS , where S is the sum of the supergravity action and the 
effective action for the D7-branes. Since A^ does not enter the supergravity action, the form of this action will not be 
needed here. Moreover, the Wess-Zumino part of the D7-branes action does not contribute for the brane orientations 
and the gauge field polarizations considered in this paper. Therefore, it suffices to consider the Dirac-Born-Infeld 
part of the D7-branes action, and since we are only interested in the two-point function in ([S]), we keep terms up to 
quadratic order in the gauge field and use the simplified form of the action 

S = —NfT D7 jf d*x ^tet~g [\ + ^§)1 F 2^ . (9) 

In this action F — dA is the overall [/(l) field strength, F 2 = g kl g mn Fk m F[ n , T D7 = 1/(2tt£ s ) 7 g s £ s is the tension of 
the D7-branes and g mn is the metric induced over them, given at zero temperature by 

ds% 7 = L- l -±f (-dxl + dx 2 ) + -^df 2 + ^dnl , (10) 

where the radial coordinate over the D7-branes is f 2 = (r 2 — L 2 )/L 2 . 
To get the equations of motion for A^ we Fourier-decompose them as 

Vz°,x,r) = J ^^ e -^ +^A M ( g ° !q ,r), (11) 

and choose q to point in the a^-direction. Under these circumstances the equations of motion for the A^-components 
are 



WM + pA )+d f (V=gg rr g 00 d f A )=o 

^g 11 g 00 w(wA 1 + P A )+d f (^g ff g 11 d f A 1 )=Q 
^-gg l \g m w 2 + 5 V)^+M^/V i 0f^)=° 



(12) 



where we have renamed q° as w, q 1 as p, and the index i takes values 2 and 3 with no summation implied over it. 

For the induce metric [TUl the first and second equations in (fT2j) imply Aq = -rj^Ai + + Ca, with C\ and C2 
integration constants. To apply the method in [l6[ and obtain H^ 11 we need the components A^ to be regular at f = 
and satisfy lim^oo A^ — 1. Given that, as we will see below, the non trivial solutions are different from zero at f = 0, 
we conclude that C\ — C2 = 0. 



4 



Using Ao = -£A\ and (JTUJ) we see that the equations p"2")) can be decoupled from each other and, furthermore, they 
reduce to the same equation, 



A + df{r J dfA) = 0, 



* (1 + f 2 ) 2 ' 

where A stands for any of the A^s and we have defined the dimensionless quantity 

~2 _ R A 2_ 9sK q 2 



Q 



L' 2 



(13) 



(14) 



The solution to (|13[) that is, as required above, regular at f — and satisfies linif-j-oo A = 1, is uniquely given by 

(15) 



a = r 



<? \ ( 3 + yrr: 



(1+0 — 5 — j? 



-;2,-f" 



2 / I 2 J \ 2 2 

where T and .F are, respectively, the Euler gamma function and the standard hypergeometric function. 



A. A brief comment on the non-zero temperature case 



In the gauge theory side we know that at zero temperature only q 2 



p 2 enters the calculation of the vacuum 



polarization function, but for the 1^0 case the solution has w and p separately as relevant parameters. Here we 
recover this behavior correctly from the gravity side, since from (|13p it is clear that the only relevant parameter 
concerning the momentum is q 2 = w 2 — p 2 , but for the T ^ case, one of the modification to the induce metric g is 
that g 0Q —g 11 , and it is not difficult to see then that the space of solutions to (jT2J) , and therefore the result for the 
vacuum polarization function, is parametrized by w and p independently. Since we are working in the T = case, we 
don't need to assume anything about the values of w or p independently, nor pick light-like momentum w — p and we 
can use q 2 as our only momentum parameter. 

B. Holographic Renormalization 

If we substitute (fTU|) and (fl~5)) into the action (|9]), the integral diverges as we extended it to f — > oo, which implies 
that we have to holographically renormalize @. It is not difficult to see that the boundary counterterm we need to 
add is 



(2ir£ 2 ) 2 , 

5S = N { Tv 7 jd /L xdQ.s —R ln(f)V-det hF 2 
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(16) 



where dD7 is the hypersurface located at r max , h a p the metric induced over it and F 2 = h afi h^ s F ai F p& . The indices 
from the beginning of the Greek alphabet run over all the coordinates of the D7-brane except for f. 



IV. COMPUTATION AND ANALYSIS 



Using © and (|15j) as an input in the method described in we obtain the expression 



iW) 



2JV ( jV c 



3tp(27r) 8 f^oo 



Jim f 6 A*(q 2 ,r)drA(q z ,r) 



(17) 



where we have already taken the trace of . 

When taking the limit in (|17|) using the actual expression (fl"5|) for A(q 2 , f), we find a logarithmically divergent term 
which, as it should be, gets canceled by the contribution of (fT5|) to the correlator in ©. The reason why we can 
explicitly check the exact cancellation of the divergences is because in the case we are analyzing, of zero temperature, 
we have analytic expressions for (fTTj) and (fl"6|). This procedure yields then a finite result for II 9 (g 2 ), that we substitute 
in ([2} to obtain 



2N t N c 



Im 



-2 + 2y / l + (q2) 



H 



1- vT + (<? 2 ) 



H 



-l + Vl + (<? 2 ) 



(18) 
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TABLE I: Comparison of the second pole prediction against experimental measurements 



Quark mass (GeV) 


Prediction (GeV) 


Measurement (GeV) 


M u = 0.266 


1.334 


1.465 


M s = 0.352 


1.765 


1.680 


M c = 1.069 


3.686 


3.771 



where H is the harmonic number function. We should notice here that H (z) has an infinite set of singular points at 
z = —k for k g N + . These points are single poles with residue -1, making R(q) a series of resonances located at 



q n = 2^n{n + 1) for n = {1, 2, 3, ...} . (19) 

We note that, although computed in a different way, the location of the resonances agrees with the spectrum of 
masses found in [l8[ for the vector mesons classified therein as the £ — case of the type II vector fields, which 
correspond to the ones we study here. The expression (fl~9]) is well defined for complex values of a 2 , so it provides the 
analytic extension to the complex plane of the distribution obtained following the method in (l8| . 

Just to illustrate the effect of the quark mass in the spectrum prediction, we obtain the quark mass by fixing the 
first pole position to its experimental value through Eq. (|19| (n = 1, N c — 3), and predict the corresponding second 
pole position. In Table |l] we compare the predicted and experimental second pole position for u, s and c -quark 
vector-meson states, using the p(770), </>(1020) and J/ip(3096) as the first poles, respectively. The prediction for the 
second pole of the u channel agrees within 9% with the experimental measurement, while for higher quark masses we 
find an agreement of around 5%. 

Given that with the method used here we obtain an analytic expression, R{q 2 = 9s ^ c ) , valid for q 2 over the 

complex plane, we can give a small but finite imaginary component e to q 2 and take its real part to be equal to s, that 
is, q 2 — s + ie. By doing this, it is possible to adjust e to account for the finite width of the mesons, and once the value 

of e has been fixed, R (q 2 — ^%(s + ie)^ becomes a function of s, R(s). This procedure is similar in spirit to the 

usage of the Breit-Wigner formula to describe unstable states in quantum field theory, suggesting that the imaginary 
part of the propagator could be seen as a product of the mass and decay width (T) of the state, namely, T = e/^/q^. 
In the large N c limit that we are working at, the mesons are fully stable and this is reflected on the fact that when 
e — > we obtain a series of delta functions , which is consistent with what has been conjectured before for this way of 
computing the spectral function [13, [H, Ho| ■ 

In order to have a qualitative insight on the predicted shape of i?(s), by including the decay width as discussed 
above, in Figure [1] we plot R(s) as obtained in the model (solid line) and the measured values (symbols), as obtained 
from e + e" scattering Q. We did choose the first pole to coincide with the p(770) vector meson mass and decay width 
parameters, and the second pole position is predicted as we discussed previously. It is worth to mention that, since 
the width is not predicted, the width exhibited for the second pole is meaningless. It corresponds to having kept the 
same width as the one used for the p(770) but substituting the predicted value for the mass of the second pole in 
the expression for the e parameter. Notice that in our model, isospin is an exact symmetry and therefore degenerate 
isospin states will not show as separate bumps. For the sake of clarity we have omitted to plot the pole associated to 
the s quark state (0). 

Experimental information on R(s) is also obtained via r hadronic decays by invoking isospin symmetry. In the 
upper panel of Figure Q] we compare R(s) (solid line) with the corresponding experimental measurements from t data 
(symbols) Q. In each case we observe a good agreement between them. In particular, the shape is very well described. 
An obvious difference between the plots is the absence of the u meson in the r data, which appears in the e + e~ data 
as a tall peak just above the p peak. This is due to the fact that the former goes through a charged current channel. 

To estimate the deviations in the shape of R(s) compared to the experimental data, we compute its contribution 
to the anomalous magnetic moment of the muon a^ l ad , 

„had ( , " 2 (0) [°° , K{s) 



3tt 2 



4m; 



This observable is suitable in this case, since the electromagnetic kernel K(s)/s [21( suppresses contributions for 
y/s > 1 GeV and therefore only the first pole is relevant. We obtain a^ ad (mr) = 660 x 10~ 10 , which is about 25% 
higher than the one obtained using experimental data (for example 520 x 10~ 10 from r data [f|). However, note that 
our value is not sharp, since the R(s) normalization is not unique. 
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0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

s 1 ' 2 (GeV) 

FIG. 1: R{s), prediction (solid line) and the measured values from e + e~ scattering (symbols). Upper panel, R(s) predicted 
(solid line) and the experimental measurement from r decay (symbols). See text for details. 

V. CONCLUSIONS 

In summary, we have provided an analytic expression for i?(s) obtained using holographic techniques at zero 
temperature in a strongly coupled, large- N c SU(N C ) field theory with N f <C N c quark flavors. We saw that our result 
for R(s) predicts the same mass spectrum for these kind of vector mesons as the one that had been established using 
other methods, and in addition it permits to include in a natural way the unstable character of the resonances. By 
comparison with the experimental data in the non-perturbative regime we can indeed recognize some of the gross 
features, as displayed in Figure [TJ in particular there is good agreement on the shape of the first pole. Since we are 
working in a theory different from QCD, the localization of the second pole is not accurate enough but improves for 
higher quark masses. In addition, the prediction of the hadronic contribution to the anomalous magnetic moment 
of the muon is well approached. It is worth pointing out that our results cannot be used to check the high energy 
behavior since the action (j9]) is dual only to the strongly coupled regime of the gauge theory. 
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